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Introduction 


The behavior of a wall jet flowing around a curved surface has been the subject of 
study for almost two hundred years following the first observations of Young* 1 * in 1800. 
Practical application of this phenomena, popularly known as the Coanda effect, have been 
pursued in more recent times with a view to delaying separation of the boundary layer on 
lifting surfaces by using a wall jet over a rounded trailing edge: the wall jet causes the 
flow to remain attached to the surface, displaces the rear stagnation point and induces 
additional circulation, and therefore produces additional lift on the surface. The use of 
this phenomena in the design of circulation control wings (CCW) has received considerable 
attention in recent years and has been described by a number of investigators (see for 
example Wood* 2 * and Englar* 3 ^). 

Despite the fact that the phenomena is being used in practical applications to wing and 
rotor design the understanding of the properties of wall jets is still limited and depends 
primarily in experimental information regarding the turbulent mixing, its effect on the 
spreading rate of the wall jet and the corresponding deceleration of the flow as it proceeds 
along the wall. A thorough review was made by Launder and Rodi* 4 ) in 1981 of available 
data for wall jets adjacent to surfaces of plane, cylinder and logarithmic spiral surfaces 
and substantial use may be made of this information to distinguish the influence of wall 
curvature, and rate of change of curvature in the stream direction, on the spreading rate 
of the jet. This information together with data from the same review on the mean velocity 
profile in the wall jet provides the basis for checking theoretical models of the flow. 

In the present paper an attempt is made to formulate the simplest possible model for 
the flow of a wall jet emanating from a two dimensional source into quiescent surroundings 
in the presence of a wall of arbitrary shape. The method uses self similar profiles for the 
mean velocity together with a simple eddy viscosity model. The streamwise and radial 
momentum equations are integrated across the wall jet flow to give an expression for the 
momentum balance including the effect of the pressure gradient induced by the rate of 
change of surface curvature. The streamwise momentum equation is also evaluated at the 
point of maximum velocity to provide a second equation and thereby permit a solution 
for the two unknown quantities b, the jet half width and u m , the maximum velocity. 
This approach provides approximate closed form solutions for the flow of the wall jet over 
surfaces of various shapes (for quiescent surroundings), and in particular permits a direct 
comparison with the available experimental results for plane, cylindrical or logarithmic 
spiral surfaces. 
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Analysis 


The wall jet is considered to comprise two parts: an inner flow adjacent to the wall 
having a highly non-linear velocity profile characteristic of a turbulent wall flow, and an 
outer flow having a velocity profile more typical of a free turbulent plane jet. 

The primary parameters that describe the flow are shown in figure 1. The jet emerges 
from a point source into a fluid at rest and spreads, increasing its width and decreasing 
its velocity due to turbulent diffusion in the jet and friction at the wall. At a distance s 
downstream of the jet exit the velocity to can be expressed as 


u = n m {a)f(y/b) (1) 

where u m is the maximum velocity, occurring at y = y m (s) , and b = 6(a) is the half 
width of the jet (at which point u = |u m ). 

The velocity profile in the outer flow (y > y m ) is assumed to take the form: 


u = u m sech 2 


(m) 


1 -& 


> & 


( 2 ) 


where ( = y/b . This velocity profile is suggested by the classical free jet solution 
found by Tollmien, modified to give u = u rn at { = £ m . The constant k is determined such 
that u = |u m at £ = l(y = 6), thus 


k = tanh" 1 = .8814 (3) 

The velocity profile for the inner flow is assumed to depend on the variable (£/ ( m ) " 
as suggested by turbulent wall flow, and is chosen to give a maximum value u = u m at 
£ = £ m : the following profile satisfies these conditions 

“ =u " [*(£)’ ”(£)’] (4) 

The value ( m , giving the location of maximum velocity, is determined by matching 
the second derivative of the velocity profiles as given by equations (2) and (4), (the first 
derivative is zero since this is the maximum velocity point). 

The result is written 


U = (l+kn)~ l 


( 5 ) 
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Before proceeding with the analysis, a comparison is made in figure 2 of the velocity 
profiles with the experimental data of Tailland* 5 )* # ). The data are typical of those taken 
by a number of authors for wall jet flows over plane and curved surfaces and show that 
the assumed profile given by equations (2) and (4) is reasonable. Additionally in figure 3 
a comparison is made of ( m , given by equation (5), with the experimentally determined 
values given by Forthmann* 7 ), Sigalia*®) Bradshaw & Gill* 9 ), Patel* 10 ) and Giles* 11 ). All the 
data falls within a band .14 < < .16 corresponding to 7 > n > 6 for Reynolds number 

R, in the range 10 4 to 10 s ; thus the values of n are consistent with those expected for 
turbulent wall flows. 

In this paper the two primary flow variables to be determined are the jet half width, 
b, and the maximum velocity Um. The continuity and momentum equations for an incom- 
pressible fluid are written 



d_ 

dy 


(hv) = 0 


( 6 ) 


. du d . . —hdp d / 2 \ 

'““57 + hv Ty {hu) = T57 + Ty ih T) 




and 


R p dy 


( 8 ) 


where 


h = 1 + y/R 

Equations (6) to (8) can be combined to give 


£.(„*) + _ h 2 uv _ L (p») \ 

ds ( } R 2 da dy\ ds \ p )> 


( 0 ) 


neglecting terms of 0 (j?)- Integration across the flow gives the integral form of the 
momentum equation: 



J_dR 
+ R 2 ds 



( 10 ) 


This equation retains the term which reflects the curvature-induced pressure gradient 
which may be large compared with r 0 /p, the wall shear stress. 
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Substitution of the velocity profiles (equations 3 and 4) into the integrals of equation 
(10) gives 


A n 




+ Cr 


^_dR 
R* ds 



(H) 


where 


An = (i +!»)(! +*»r l 
= \ [(| lo « 2 - 1 ) ” 2 + 1 -+ l ] (> + 

FVom equation (11) it can be seen that in the absence of curvature effects the wall jet 
momentum {bu 2 m ) decreases with distance due to skin friction at the wall. However for a 
wall of decreasing radius of curvature (" < 0) it is possible for the wall jet momentum 
to increase with distance along the wall. This will occur if 



> 



An approximate form of equation (11) valid for large n and for C/-*0 may be written 


- * ds K m) k R 2 dS 


U: 


(llo) 


a result which is applicable for large Reynolds number. 

Since there are two unknown quantities in equation (11) it is necessary to use a second 
relationship in order to determine b and tim. It is convenient to use the momentum equation 
(7) evaluated at the point of maximum velocity, (ie where u — u m , §* = 0) which may be 

written 


du uv 1 dp 


hdr 2 r 

U 3a + ~R + pds 

m 

p dy + pR 


(12) 


The terms on the left hand side may be written in dimensionless form, neglecting 
terms of 0{b/R ) 2 as follows: 


102 


and 



Finally, 


6 1 dp 

u 2 m ' pda 



|n b \ 1 d fL 2 . bdR 
1 + kn ' R '[u 2 m d3^ bUm ^ Rds 

Thus the left hand side of equation (12) becomes 


1 

2 



1 + fn 
l + *n 




1 

1 + kn 



2 / n \ b 2 dR 
+ 3 Vl + kn) R 2 da 


In the limiting case of large n this expression simplifies to 



db 4_ 
da + 3k 



In order to evaluate the right hand side of equation (12) it is necessary to express r in 
terms of u through an eddy viscosity e based on a length scale 6(1 — £ m ) and a characteristic 
velocity u m . It is convenient to write 


6 


£fill - fm)» «. ,(() 


(13) 
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where K is a constant to be determined experimentally. The function g{£) is chosen such 
that y(6») = 1 and also satisfies the appropriate boundary condition at £ = O. 

The shear stress takes the form 


T _ fdu tt \ 

~p- € [dy R)’ 

or, substituting for « = u m /(0 and e 



£(i -U)\f9-b/Rf9] 


(M) 


(15) 


Similarly, 


b dr 

pu 2 m dy 


- «-) [(/'» + /V) - 1(/ , » + N) 


(16) 


where / denotes a derivative with respect to £. 

Since f{£ m ) = g{U) = 1 and f'{£ m ) = O the foregoing expressions evaluated at 
£ = £ m combine to give: 


(-~ + ~) = T 1 { J?(l - f»)(! + f- ‘/*)£(1 - »/* s'//"). + O(V-R) 1 } (17) 

pf2/ m 6 I4fc 2 J 

The quantity /" is evaluated from 


/ = sech 2 


~ £»») 

. . 


as 


ru«.) = 


2* 2 

(1 - *m) a 


Turning now to the function y(£) which describes the variation with £ of the eddy 
viscosity, y(£) is assumed to take the form: 




( 18 ) 
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that 


where the exponent m and the constant a are determined by satisfying the condition 


pu; 


Substitution for c and u gives 


t (du u\_l 

•4 e U 2 C,a( 


K 1-6. „ (±\ 

(. 1 UJ 


m+i-1 




at £ = o. 
Thus 


, 1 , C/* 

m = 1 and a = — £— 

n K 


(since ^4® = kn) 

Thus 

and tfitm) is evaluated as 



, w 

|JL' 

UJ 1*1 

► 

1 

>! 

6». 


= (2 - i - (1 + *») 


Thus equation 17 becomes 


(19) 


( 20 ) 


(20a) 


2b (hdr 2r\ kn f. . [* “ n / 2 + (l - fir)] 

<{-pr y + jR) m — tth { 1+ TTkH hlR \ 


which has the form, for n large and Cj — ► 0: 


26 (hdT 2T\ _r n . ._1 


( 21 ) 
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Substitution into equation (12), and elimination of the quantity ^r£(6t4) using equa- 
tion (11a) gives finally, an expression for the spreading rate g: 

£=4 + -M + (If - - ¥ ) ( i/fl f ) b/R + ° WR) ' (22) 

In this equation K is to be determined experimentally for the flat plate (R -*• oo) 
and the term involving (b/R%) is retained since this may be the same order as g (for a 

log spiral shape, in fact (fc/ii^) = 

It is clear from equation 22 that the influence of small curvature (small b/R) may be 
quite significant since the quantity £ • b/R appears in the expression for the spreading 
rate. 


Comparison with experimental Results 

In order to validate the foregoing analysis a comparison is made of the spreading rate 
as given by equation (22) with available experimental results for two cases, namely, the 

circular cylinder (R = constant) and the logarithmic spiral = const ). 


(a) Circular Cylinder 

In this case equation 22 reduces to a linear differential equation for b/R namely 

(22a » 

The solution to this equation may be written in the form 


» K i-v* 

« log(l + fV«) 


or, neglecting terms of 0(b/R) 7 , 


(23) 


6/s = K (l + ^ b/R) (23a) 

The constant K is determined from the flat plate wall jet (b/R = O) as K = .07, and 
equation (23a) with k = .8814 gives 
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6/s = 07(1 + 46/12) 


(n = 7) 


(236) 


6/5 = .07(1 + 3.46/J2) (n = 6) (23c) 

These expressions are shown in figure 4 together with the available experimental re- 
sults of Fekete^ 12 ), and Wilson and Goldstein* 13 *. In this regard the review by Launder and 
Rodi* 14 * suggest that the data of Fekete is possibly more representative of two dimensional 
flow than that of Wilson and Goldstein. It is in fairly good agreement with the theoretical 
expression for n = 6 (the Reynolds number is in the range 4 to 13 x 10 s ). 

(b) Logarithmic Spiral 

For a logarithmic spiral R/s = constant and equation (22) has a self similar solution 


6/s = 


[l _ Itplt/R 


rj 07 |1+86/Jl1 

* 07 ii- 6 7*r 

.07[1 + 6.86/R 

* [1 — 6 /R] 


(n = 7 
(n = 


(24) 


(24a) 

(246) 


A comparison of these curves with the experimental results of Guitton and Newman* 14 * 
and Kamemoto* 15 * is made in fig. 5. Here Guitton and Newman’s results are considered to 
be most representative of two dimensional flow. Equation (24) reflects the much greater 
increase of b/s with b/R for the logarithmic spiral, compared to the circular cylinder 
equation (23a), and this trend is also seen in the experimental data. 

Conclusions 

FVom the foregoing analysis and comparison with experimental results it is seen that 
the wall jet is influenced by three effects: 

(a) The turbulent diffusion in the wall jet giving rise to a linear spreading rate db/ds = 
K « .07. 

(b) A coupling between the eddy viscosity and the curvature of the wall arising from the 
term J^( eu/R ) in the stress gradient and yielding a contribution of £ = Kjb/R. 
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(typically this contribution is 0{\b/R) for Reynolds numbers typical of laboratory 
experiments), and 

(c) A curvature-induced streamwise pressure gradient giving rise to a contribution of 
0(b/R) 2 %. When g < O this effect causes the jet momentum (6t&) to increase in 
the streamwise direction and help confine the jet to a thin layer as it proceeds along 
the surface. This is an inertial term, largely independent of Reynolds number, and is 
an essential feature of Coanda turning of the flow. 

These influences are additive and give an approximate relationship between the spread- 
ing rate g, the half width of the jet and the surface radius of curvature R, of the form 





dR 

da 


which can be integrated to give b(s) when the wall shape R(s) is known. 

The use of this result for the circular cylinder (R = const) shows good agreement 
with the previous experimental results of several authors. This provides some confidence 
that a simple algebraic eddy viscosity model, used in conjunction with the appropriate 
expression for the shear stress (including curvature terms), is sufficient to describe the 
primary features of the wall jet. 
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Figure 2. Comparison of assumed velocity profile with experimental data. 
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